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Abstract 



The question of the optimality of Thompson Sampling for solving the stochastic multi-armed 
bandit problem had been open since 1933. In this paper we answer it positively for the case of 
Bernoulli rewards by providing the first finite-time analysis that matches the asymptotic rate given in 
the Lai and Robbins lower bound for the cumulative regret. The proof is accompanied by a numerical 
comparison with other optimal policies, experiments that have been lacking in the literature until 
now for the Bernoulli case. 

1 Introduction 

In a stochastic bandit problem an agent is repeatedly asked to choose one action from an action set, each 
of which produces a reward drawn from an underlying, fixed, but unknown distribution associated with 
each action. Thus he must choose at each time whether to use the observations he has already gathered 
to gain the greatest immediate reward (exploitation) or whether to choose an action from which few 
observations have been made and risk immediate loss for greater knowledge and potential future gain 
(exploration). In this paper we focus on stochastic bandits with Bernoulli rewards, initially proposed 
by Thompson in his paper of 1933 [14] to model medical allocation problems. Thompson's paper also 
presented the first bandit algorithm, Thompson Sampling. This algorithm has received much attention 
in the recent literature, and in this paper we give the first theoretical proof of the asymptotic optimality 
of this algorithm in the context of cumulative regret minimisation. Furthermore we achieve this result 
by giving a finite time analysis for the algorithm. 

Associated with each action, a, is an unknown Bernoulli distribution B (fJ- a ), whose expectation is fi a . 
At each time t the agent chooses to observe an action A t <E {1, . . . , K} and receives a reward Rt drawn 
from the distribution B (fJ-A t ) ■ A policy, or bandit algorithm, is defined to be a (possibly randomised) 
method for choosing A t given the past history of observations and actions. The agent's goal is to minimize 
the expected cumulative regret of his policy, which is defined to be: 



where fi* = max a /i a denotes the expectation of the best arnj^j or optimal action, and N a< t the number of 

draws of arm a at the end of round t. Lai and Robbins proved in |10) that all strongly consistent policies 

x The words arms and actions are used interchangably. 




(1) 




(i.e. policies satisfying lZ(t) — o(t a ) for all a G (0, 1)) must satisfy, for any suboptimal arm a 

liminfg^> , 1 , (2) 

where K(p,q) denotes the Kullback-Leibler divergence between B(p) and B (q): 

V 1 — V 

K(p,q) :=pln- + (1 -p)ln- . 

q l-q 

Their result, which holds for more general classes of reward distributions, leads to the definition of asymp- 
totically optimal policies as policies that satisfy ^ with equality. 

In the same paper [lOj Lai and Robbins were able to describe an asymptotically optimal policy, 
however no finite-time analysis was provided, nor was it an efficient policy to implement. The UCB1 
algorithm by Auer et al. [4] was the first of a series of efficient policies, like UCB-V [3] or MOSS [2], for 
which good regret bounds in finite time were also provided. These policies all use an upper confidence 
bound for the empirical mean of past rewards as an optimistic index for each arm, choosing at each 
time the action with the highest current index. However, for each of these algorithms we only have the 
result that there exists two constants K\ > 2 and K > such that for every suboptimal action a, with 
A a = /i* - |Lt 0) 

E[7V a ,T] < ^\n(T)+K 2 . (3) 

This does not imply ^ with equality since by the Pinsker inequality 2K(fi a , fj,*) > A^. On the contrary, 
recently proposed index policies such as DMED [8] and KL-UCB [6l [IT] , which use indices obtained from 
KL-based confidence regions, have been shown to be asymptotically optimal. 

Unlike most of this family of upper confidence bound algorithms that has been so successful, Thomp- 
son Sampling is a policy that uses ideas from Bayesian modelling and yet it solves the fundamentally 
frequentist problem of regret minimisation. Assume a uniform prior on each parameter fj, a , let 7r ajt denote 



the posterior distribution for /i a after the t round of the algorithm. Let 9 a j denote a sample from ir a .u 
we sometimes refer to 9 a t as a Thompson sample. Thompson sampling is the policy which at time t 
chooses to observe the action with the highest Thompson sample 8 a t , i.e. it chooses action a with the 
probability that this action has the highest expected reward under the posterior distribution. 

Before Agrawal and Goyal's recent paper pQ Thompson Sampling had been investigated in [7] as the 
Bayesian Learning Automaton, and in |12j where an optimistic version was also proposed; however these 
papers only provided weak theoretical guarantees. In [3] extensive numerical experiments were carried 
out for Thompson Sampling beyond the scope of the Bernoulli bandit setting (to the Generalized Linear 
Bandit Model) but without any theoretical guarantee at all. Consequently the first finite-time analysis of 
Thompson Sampling in [TJ was a major breakthrough, yet the upper bound for the regret that is shown 
in this paper scales like ^ and the question of Thompson Sampling's asymptotic optimality was still open. 

Meanwhile, there has been a resurgence of interest in Bayesian strategies for bandit problems (see [H] 
for a review of them). The Bayes-UCB algorithm, an upper confidence bound policy which uses an adap- 
tive quantile of ir a ,t as an optimistic index, was the first Bayesian algorithm to be proved asymptotically 
optimal. In this paper we are able to show that the same is true for a randomised Bayesian algorithm, 
Thompson Sampling. Moreover we refer in our analysis to the Bayes-UCB index when introducing the 
deviation between a Thompson Sample and the corresponding posterior quantile. 

Contributions We provide a finite-time regret bound for Thompson Sampling, that follows from ([I]) 
and from the result on the expected number of suboptimal draws stated in Theorem [2j More precisely 
we show the following: 
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Theorem 1 For every e > there exists a problem- dependent constant C(e, fj,i, . . . , /iR-) such that the 
regret of Thompson Sampling satisfies: 



K(T)<(l + £ ) £ A.(ln(r) +h ,l„(T)) +cfe(ii w) _ 



Besides this asymptotically optimal regret bound, we also provide the first numerical experiments that 
show Thompson Sampling outperforming the current best optimal policies like DMED, KL-UCB or 
Bayes-UCB. The rest of the paper is structured as follows. Section [2] contains notations or results already 
used in PP, [B] or [5j that are useful in our finite-time analysis given in Section [3] Numerical experiments 
are presented in Section |4j 

2 Preliminaries 

We gather together here some useful preliminaries such as notations not already given in the introduction: 

• For the rest of this paper, we assume action 1 is the unique optimal action. Without loss of 
generality we can assume that the parameter /i = (/xi, of the problem is such that [i\ > 
A*2 > ••• > li-K- 

• We shall denote by S a ,t the number of successes observed from action a by time i, and denote the 
empirical mean by: 

_ S a j 

• In the Bernoulli case, with a uniform prior on the parameters /i a of the arms, the posterior on arm 
a at time t is explicitly 

ir a ,t = Beta (S a , t + 1, N a , t - S a ,t + 1) • 

• Let g ta denote the cdf of a Beta distribution with parameters a and b and (resp /? ) the 
cdf (resp pdf) of a Binomial distribution with parameters j and fi. We recall an important link 
between Beta and Binomial distribution which was used in both [T] and [5]: 

F»t(y) = l-F* +b _ 1 Ja-l) 

We use this 'Beta-Binomial trick' at several stages of our analysis. 

• We denote by u a ,t (resp. q a .t) the KL-UCB (resp. Bayes-UCB) index at time t, and define them by 

(rsfSa.t \ . ln(t) + ln(ln(T)) | 

X> JV a ,( 

where Q(a, n) denotes the quantile of order a of the distribution tt 

• A special link between these two indices is shown in [S]: q a t < u a t . 

2 In Appendix A of pQ the authors show that adding a second optimal arm can only improve the regret performance of 
Thompson Sampling. 
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3 Finite Time Analysis 



3.1 Sketch of Analysis 

Unlike Agrawal and Goyal's analysis, which is based on explicit computation of the expectation EfiVjjr], 
we are more inspired by standard analyses of frequentist index policies. Such policies compute, for 
each arm a at round t, an index, l at , from the sequence of rewards observed from a by time t, and 
choose A t = argmax a Z at . A standard analysis of such a policy aims to bound the number of draws of a 
suboptimal arm, a, by considering two possible events that might lead to a play of this arm: 

• the optimal arm (arm 1) is under-estimated, i.e. < Mi! 

• the optimal arm is not under-estimated and the suboptimal arm a is drawn. 

Taking these to be a good description of when the suboptimal arm is drawn leads to the decomposition 

T T 

E[N a , T ] < ^ P (l ht < m ) + P ((*«*,* > h,t > Ml) n {A t = a)) 
t=i t=i 

The analysis of an optimistic algorithm then proceeds by showing that the left term (the "under- 
estimation" term) is o(ln(T)) and the right term is of the form K , > ln(T) + o(ln(T)) (or at worst 

^2 ln(T) + o(ln(T)) as in the analysis of UCB1). This style of argument works for example for the 
KL-UCB algorithm [6] and also for the Bayesian optimistic algorithm Bayes-UCB [S]. 

However we cannot directly apply this approach to analyse Thompson Sampling, since the sample 
9 a ,t is not an optimistic estimate of \i a based on an upper confidence bound. Indeed, even when ti\^ is 
well concentrated and therefore close to a Gaussian distribution centred at Hit P (#x,t < Mi) is close to | 
and the under-estimation term is not o (ln(T)). Hence we will not compare in our proof the sample 9 a t 
to fi a , but to n a — \n(t)/N a j (if jV 0ii > 0) which is the lower bound of an UCB interval. We set the 
convention that if N a j. = 0, y/61n(t)/N a j — oo. 

As is observed in £Q the main difficulty in a regret analysis for Thompson Sampling is to control the 
number of draws of the optimal arm. We provide this control in the form of Proposition [T] whose proof, 
given in Section [3.3| explores in depth the randomised nature of Thompson Sampling. 



Proposition 1 There exists constants b = b([ii,[i2) € (0, 1) and Cb = Cb()J-i, H2) < 00 such that 

00 

t=l 

This proposition tells us that the probability that the algorithm has seen only a small number draws 
on arm 1 is itself small. As a result we can reduce to analysing the behaviour of the algorithm once it has 
seen a reasonable number of draws on arm 1, and thus the posterior distribution is well concentrated. 

Remark 1 In general, a result on the regret like E[Ni, t ] > t—Kln(t) does not imply a deviation inequality 
for N\ t i (see }13f )- Proposition^is therefore a strong result, that enables us to adapt the standard analysis 
mentioned above. 

Using this result, the new decomposition finally yields the following theorem: 

Theorem 2 Let e > 0. With b as in Proposition^ for every suboptimal arm a, there exist constants 
D(e, Hi, fi a ), N(b, e, Hi, [i a ) and N (b) such that: 

E[JV 0lT ] < (1 + e ) ln ^) + lnln ( T ) + D ( e Mi, Ma) + N(b, 6, mi, Ma) + N (b) + 5 + 2C b . 
iv(Ma,Ml) 

The constants will be made more explicit in the proofs of Proposition [T] and Theorem [2] The fact that 
Theorem [2] holds for every e > gives us the asymptotic optimality of Thompson Sampling. 
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3.2 Proof of Theorem [2] 

Step 1: Decomposition First we recall the modified decomposition mentioned above: 



<^P (^, t < Ml 



'61n(t) 
N ljt 



The sample # Qjt is not very likely to exceed the quantile of the posterior distribution q a j we introduced: 

|:^..<>^<|^<i^p<2 

where this last inequality follows for T > e. So finally, using that u at t > q a .t. 



E[N a , t ] <j2p U t < mi - y^^l + E P ( M <^ > w 
"> ' "> 



'61n(t) \ „ 

W -,A t = a +2 (4) 



i.i 



Step 2: Bounding term A Dealing with term A boils down to showing a new self-normalized in- 
equality adapted to the randomisation present in each round of the Thompson algorithm. 

Lemma 1 There exists some deterministic constant No(b) such that 



E P rM<^i- t/^rr J <N (b) + 3 + C b <oo 



with Cb defined as in Proposition^ 



Proof. Let (Ut) denote a sequence of i.i.d. uniform random variables, and let Ei >s be the sum of the 
first s rewards from arm 1. In the following, we make the first use of the link between Beta and Binomial 
distributions: 

U t <l-F* ( ^ QS M ) I n (N ht > t b ) I + P (N ht < t v 



F U r— - (5i jt ) <U t \n (N lit >t b )\+¥ (N ht < t b ) 



Afi.t+l.Mi- 

< p (as e {£"...;} : ^ B +i / ^ (£ M ) < ^ + p (JVi, t < * b ) 

= E P ( S M < J + P {N ut < t b ) 



■5 



The first term in the final line of this display now deals only with Binomial random variables with large 
numbers of trials (greater than t b ), and so we can draw on standard concentration techniques to bound 
this term. Proposition [T] takes care of the second term. 



Note that (F B ) 1 , lUt) ~ Bin ( s + l,Ui — \/6 \n(t)/s ) and is independent from Ei 

s+l,[i 1 -y/61n(t)/s \ J 

Bin (s, fii). For each s, we define two i.i.d. sequences of Bernoulli random variables: 



(*!,,)! ~ B I Mi - \ I and (X a ,,)i ~ 6 ( m ) , 



and we let := X 2 .j — -^u, another i.i.d. sequence, with mean \/ 6 ■ Using these notations, 



>/*l— V 



vi=i / \;=i 



Let iV Q (&) be such that if f > AT (&), V 6 ^ M*) ~ 1 > V 5 ^ ln M- For 1 ^ ^o(&), we can apply Hoeffding's 
inequality to the bounded martingale difference sequence Z[ = Z% — W 6 ln(t) / s to get 



We conclude that 



Step 3: Bounding Term B We specifically show that: 

Lemma 2 For all a = 2, . . . , K , for any e > there exist N(b, e, fix, fj, a ), L>(e, /Xi,/x a ) > suc/i i/iai /or 
aH T > jV(6, e, /ii, /i a ) 

, N , ln(T) + lnln(T) 

KyfJ-a, Ml) 

Proof. First rewrite term J5 so that we can apply Proposition [T] 



t=l 

T 



t=i \ 



b 



For ease of notation we introduce 

K+(x,y) := K(x,y)l (x < y) , f T (t) := ln(t) + ln(ln(T)) 

/61n(t) Jr ^ . . . ln(T) +lnln(T) 

V i 6 K(^ ,Mi) 

Now 

(«o,t > a) = (iV 2 , t iir + (M2,JV 2 , t ,a) < /r(*)) 



G 



and so summing over the values of A^t and inverting the sums we get 



^PK,t >Mi- 0t,A t = a) =E 



YK T ,a\ T 

E Ei 



( S /f+(Aa, a ,Ml-ft)</T(t)) 1 (^t=a,Af 2it = S ) 



s— 1 i — s 

T T 



•E 



E E 1 (sif+(Aa, s ,A'i-/3 t )</T(t)) 1 (A t =a,Ar 2 , t = s ) 

S=L-fT,aJ+l t = S 



As y n- _K" + (/t a .s, y) is increasing and 1 1— > j3 t is decreasing for i > e 1 / 6 , for T such that 

#T,a(e) > e 1 /", 

we have that if t > Kx, a (€), 

1 (sK+(Aa,»,Mi-A)</T(*)) ^ 1 ( S K+(/ia, s ,Mi-/3K T _ a )</r(t)) - 1 (s-ff+(/ia, s ,Mi-fe T _ a )</r(T))- 



(5) 



and therefore, 

T 



Y,V(u a ,t>»i-Pt,At = a)<E 



t=i 



E E 1 (A t =a,N 2 , t =s) 



-E 



S—l t—S 

T 

E i 



(sK+(jl a . s ^ 1 -l3 KT a )<f T {T)) 

S=\K T , a \+l t=S 



E 1 (A t =a,N 2:t = 



Given that X)t=s l(At=a,JV 2 t=s) < 1 f° r au s i the nrs t term is upper bounded by K a ,Ti whereas the second 
is upper bounded by 



E 



E 1 (^T,a-R"+(Aa, S ^l-/3A- Ta )</T(T)) 



So, for T satisfying |5]), 



{B)<K T , a + J2 P U + (Aa, s , Mi - Pkt.J) < Ml) 

L-f^T.aJ+1 ^ 



Using the convexity of if + (fi a .s, ■), we can show that 



K + ((i atS ,fj,x) < AT + (/i QiS ,/ii - (3 KaT ) H v. \$K aT - 

MiU-MiJ 



If K + (Aa, s ,Mi - fex.J < K (Va, Mi)/! 1 + e), tncn 



K + (fla,s, Ml) < "Y'""" 1 ' + 7T^ T 

1 + e /ii(l-^i) 



< 



l + e/2 



(0) 



where the last inequality ^ holds for large enough T. There exists a deterministic constant N = 
N(b, e, m, n a ) such that for all T > N both ^ and ([6} are satisfied. Hence, for all T > N 

(B)<K T>a + £ p(V( AaiSlMl )<feMY 

L^T,„j+l V i+ 2 / 

Since this last sum is bounded above explicitly by some constant D(e, /xi, /x ) in [TT] we have proved the 
lemma. To be explicit, D(e, Mi, Ma) — ( 1+£ / 2 ) ■ 



2 (min(/j a (l-/i a );^ii(l-/Ji))) 



7 



Conclusion: The result now follows from Lemmas [I] [5] and inequality 
3.3 Proof of Proposition [l] 

Since we focus on the number of draws of the optimal arm, let Tj be the occurrence of the j play 
of the optimal arm (with r := 0). Let £j := (t j+1 — 1) — Tf this random variable measures the num- 
ber of time steps between the j th and the (j + l) th play of the optimal arm, and so ^2^ =2 N a ,t = Y^f=o £r 

For each suboptimal arm, a relevant quantity is C a = (jttz^ vs and let C — max a ^! C a — ^zr^yz- 
We also introduce S a = Atl ~ /J " and let 6 = 62- 

Step 1: Initial Decomposition of Summands First we use a union bound on the summands to 
extract the tails of the random variables £j : 



P(JVi,i < t b ) = P VJ N a ,t > i 




<P(3je{0,..,^J}:O>^"-l) 

L* b J 

<Y, P ^>t 1 - b -l) (7) 

This means that there exists a time range of length t 1 — 1 during which only suboptimal arms are 
played. In the case of two arms this implies that the (unique) suboptimal arm is played \- — ^] ti mes 
during the first half of this time range. Thus its posterior becomes well concentrated around its mean 
with high probability, and we can use this fact to show that the probability the suboptimal action is 
chosen a further \- — ^] times in a row is very small. 

In order to generalise this approach we introduce a notion of a saturated, suboptimal action: 

Definition 1 Let t be fixed. For any a 7^ 1, an action a is said to be saturated at time s if it has been 
chosen at least C a m(i) times. That is N a>s > C a ln(i). We shall say that it is unsaturated otherwise. 
Furthermore at any time we call a choice of an unsaturated, suboptimal action an interruption. 

We want to study the event Ej = {£j > t 1_b — 1}. We introduce the interval Tj = {-n,-, Tj + — 1] } 
(included in {rj,Tj + i} on Ej) and begin by decomposing it into K subintervals: 







'1(^-1)' 


} 


K 


,Tj + 


K 





I = 1,...,K. 



Now for each interval 2^ j, we introduce: 

• Fj f. the event that by the end of the interval Ijj at least / suboptimal actions are saturated; 

• rij i. the number of interruptions during this interval. 

We use the following decomposition to bound the probability of the event Ej : 

F(Ej) = F(Ej n Fj t K-i) + F(Ej n F!f K _ 1 ) (8) 

To bound both probabilities, we will need the fact, stated in Lemma [3j that the probability of 9i s 
being smaller than ^2 + 6 during a long subinterval of Ij is small. This follows from the fact that the 
posterior on the optimal arm is always Beta(5i. Tj + l,j — Si tT . + 1) on Xj\ hence, when conditioned on 
Si. Tj , ^i,s is an i.i.d. sequence with non-zero support above /12 + 5, and thus is unlikely to remain below 
fi2 + S for a long time period. This idea is also an important tool in the analysis of Thompson Sampling 
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Lemma 3 3Ao = Ao(miiM2) > 1 such that for A g]1, Ao[, for every (random) interval J included inXj 
and for every positive function f , one has 

P ({Vs G J, 9 1>s <(i2 + 5}n {\J\ > f(t)}) < (a MliMa )'<*> + C Ao , A yi x e-^.«-« 
where C\ imifl3 ,d\ lfllttl , > 0, and a Ml)|t2 = (1/2) 1- ' 12- ' 5 . 

The proof of this important lemma will be postponed to section |3.4| and all the constants are explicitly 
defined there. Another keypoint in the proof is the fact that a sample from a saturated suboptimal arm 
cannot fall too far from its true mean. The following lemma is easily adapted from Lemma 2 in pQ. 

Lemma 4 

P (3s < t, 3a 1 : 6 a . s > ^ + 5 a , N a , s > C a ln(t)) < 2 ' A " ' ' 



K 



we 



i 2 ' 

Step 2: Bounding P(-Ej PI -Fj,fsr— l) On the event iy n F^k-X: only saturated suboptimal arms are 
drawn on the interval Ij,K- Using the concentration results for samples of these arms in Lemma [4j we 
get 

¥{Ej n Fjvk-i) <P({3s G /,,./,.« ^ 1 : a , s > Ma + 5} n £y n F i)JC -i) 
+ P({Vs G Xj-k, a ^ 1 : a ,s < Ma + <U H £y fl 
<P(3s <t,a^l: 6> a , s > ^ a + 5 OJ iV 0)t > C Q ln(t)) 
+ P({Vs G Xj, K , a^l: 6 a , s < ^ 2 + 5} n E s n fj.jr-i) 

< 2( ^ 2 ~ 1} < i J L2 + S,Wsel j , K ). 

The last inequality comes from the fact that if arm 1 is not drawn, the sample 6\ tS must be smaller than 
some sample a s and therefore smaller than M2 + S. Since Ij k is an interval in X, of size 
get using Lemma |3j for some fixed A G]l, Aq[, 

h, s < M2 + S,Vs el jlK }) < (<wJ K +C A)Ml)P3 l — — — J e ^.w.^ = : g^,^, b, j, t). 

(9) 

Hence we have show that 

v(E 3 nF,, K ^) < 2( ^2 +g(jn,i*2,b,j,t), (10) 

and choosing b such that 6 < 1 — i, the following hypothesis on g holds: 

t>i j^t 6 

Step 3: Bounding P(-E7j fl iy c #;— l) We show through an induction that for all 2 < I < K, if t is 
larger than some deterministic constant iV^ 1)jU2) i, specified in the base case, 

P(£, n Fjj_ x ) < (I - 2) ( 2( ^~ 1} + /(Ml, M2, 6, J, t) 

for some function / such that Xit>i Si<j<t b /(Mi> A*2j b, j, t) < oo. For I = K we get 

P(^ n Ff^.i) < (K- 2) ( 2( ^~ 1} + /(mi,M2,6, j,*)) • (11) 
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Step 4: The Base Case of the induction Note that on the event Ej only suboptimal arms are 
played during Ij t i- Hence at least one suboptimal arm must be played |" - g 2~ 1 ] times. 

There exists some deterministic constant -/V MliA12 j, such that for t > A^ liM2j b, |~ * ^2~ 1 ] > Chx(t) (the 
constant depends only on \i\ and /i2 because C = C2). So when t > N^^fi, at least one suboptimal 
arm must be saturated by the end of 2j,i- Hence, for t > N^^fi 

¥{E 3 ni^) = 0. 

This concludes the base case. 



Step 5: The Induction As an inductive hypothesis we assume that for some 2 < / < K — 1 if 
* > N^^b then 

p(f, n F?,^) < (i - 2) ^ K ~ l) + f(vi, f i2,b,j,t) 
Then, making use of the inductive hypothesis, 

F(E j n Fjj) < F(E j n Fl^) + p(Ei n i^j n F,-,j_i) 

< (/ - 2) ( 2( ^ 2 ~ 1} + /(Mi.A*a,6.i,*)) + n F?, n 

To complete the induction we therefore need to show that: 

F(Ej n F?, n Fj^i) < 2{K t ~ l) + f(^^ 2 ,b,j,t). (12) 

On the event (-Ej n F, c ; f~l Fj t j_i), there are exactly Z — 1 saturated arms at the beginning of interval 
and no new arm is saturated during this interval. As a result there cannot be more than KCln(t) 
interruptions during this interval, and so we have 

P(F, n Ff, n < F(Ej n F jtl _ t n {n,- ; < KCln(t)}). 

Let Si denote the set of saturated arms at the end of X^i and introduce the following decomposition: 

P{E J n Fj t i_i n {n jt i < KC\n{t)}) 

<F({3s e %i,o e : 6» a , s > /i a + 5 a } n Ej n 

s ^ 

+ p({Vs e Jj.1,0 e : 6» a , s < Ma + <U n ^ n F 3 -,i_i n {n w < ^cin(t)}) . 

S v ' 

B 

Clearly, using Lemma [4j 

2(K-1) 



(A) <F(3s<t,3a^l: 6 a , s > Ma + S a , N as > C a ln(i)) < 



t 2 



To deal with term (B), we introduce for k in {0, . . . , rijj — 1} the random intervals Jk as the time range 
between the k th and (k + l) st interruption in Ijj. For k > rijj we set Jk = 0- Note that on the 



event in the probability (B) there is a subinterval of Tjj of length J K2 during which there are 

no interruptions. Moreover on this subinterval of Ij t i, for all a 7^ 1, 6 a s < /12 + #2- (This holds for 
unsaturated arms as well as for saturated arms since their samples are smaller than the maximum sample 
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of a saturated arm.) Therefore, 



t 1 - 6 - 1 



(B)<F[{3ke{0,...,n jtl }:\J k \> CK2Ht) 



KC\n(t) 

± E ' 

k=l 
KC\n(t) 

* E ' 

k-1 



\Jk\ > 
\Jk\ > 



1-6 



1 



CK 2 ln(t) J 

t l - b -l \ 
CK 2 ln(t) J 



n {Vs e 2j- /; a e S^l : 6> a , s < ^ 2 + 8} n £y n F^-i 
n{Vse J t ,a^l: 6» a , s < ^ 2 + 6} n J5j 

□ {VsGjfc, M <M2+<5}) (13) 

1 included 



Now, we have to bound the probability that 9\^ s < to + <5 for all s in an interval of size CK2 ln ( f ) 
in Ij. So we apply Lemma [3] to get: 



1— b_ 
^1 



(B) < CK\n{t) ( a(tllMj )^W + C A , M1 



\CK 2 ln(t; 



-i«ix,Mi.M := f( flljl j, 2 ,b,j,t). 



Choosing the same b as in (l9|), we get that Et>i Ei<j<t l f{^ii^2ib,j,t) < +oo. It follows that for this 
value of b, (12 1 holds and the induction is complete. 

Step 8: Conclusion Let b be the constant chosen in Step 2. From the decomposition (JsJ) and the two 
upper bounds (10) and (11), we get, for t > N^^^y. 



¥(E J )<(K-2) 



2(K-1) 



f(fj,i,fi 2 ,b,j,t)) 



2(K-1) 



fJ.2,b 



Recalling (|7|) , summing over the possible values of j and t we obtain: 

t b 

$>(^M < t") < N^ 2 , b + 2{K - If J2 jkb + E E^/0^ ^2, b, j,t) + gfa, to, b,j, t)] < C Ml 

t>l t>l t>\ j = l 

for some constant C Mlj ^ 2! b < oo. 
3.4 Proof of Lemma |3] 

On the interval the J (included in Xj by hypothesis), the posterior distribution 7Ti )S = 7Ti !T . is fixed and 
0i, s is, when conditioned on Si tTj , an i.i.d. sequence with common distribution Beta(5i, T . +1, j — Si jTj +1). 
Hence, 

P (6i, s <H2+6\sEj) = F^ +lj _ Si Tj+1) (to + S) = 1- F§ +1:fl2+s) (S hTi ) 

where we use the link between the tail of Beta and Bernoulli distribution mentioned above. Using the 
independence of the #i, s gives 



P (V, G J, 6 hs <to + 6\S hT] ) = (l - F§ +1 ^ +s) (S hTj )) lJl < (l - F { B ]+1 ^ +5) (S hT] )) 
Finally 

p (Vs e J, 9i,s < m + S) = e [p (Vs e J, 9i,s <to + s\s hTj )] < E (l - r, 

An exact computation of this expectation leads to 



/(*) 



(5l,r 3 ) 



/(*) 



E 



s=0 
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To simplify notation, from now on let y = f/,2 + S. Using, as in PQ, that FE_ X (s) = (1 — y)F?(s) + 
yF? y (s-l)>(l-y)F* y (s),weget: 

i 1 F (f+i.,)( s )) /(t) ^ ex p (-/W F (f+i^)( s )) ^ cx p (-/WC 1 - 



Therefore, 



E 



(1 - if 



s=0 



Using the fact that for s > \yf\ , F^ y (s) > | (since the median of a binomial distribution with parameters 
j and y is [j/j] or [yj\), we get 

E[(l-^ + i„ 2+ ,)(^)) /(i) 



by J 



<^exp(-/(t)(l- y )F ( f >2/) ( S ))/^ 1 ( S )+ £ (J 
< £ exp (-/(*)(! - y)F { l y) (s)) ff^ (s) + - 



(i-y)/W 



s=0 



(l-j-i/m 

2 



It is easy to show that for every A > l,Vx > 0, x x exp(— x) < (^) This allows us to upper-bound the 
exponential for all A > 1, using C\ — (^) A ,by: 

(E)< 



< 



(/(*)d -y)r s ( %)(a) )- - (/(t)d-y)r s (/^ )(S ))- 

Now, inspired by Agrawal and Goyal's work (proof of Lemma 3) we compute: 



e s )M?(i-Mi) j - 



"J-Fl ' \sJ ' J- v r - L/ < r-ir r~+y 



(i-y)V Vy A (i-Mi) 

Let R\(u,±,y) = ^tttt \ ■ There exists some Ai > 1 such that, if A < Ai, R\ > 1. More precisely, 



R\ > 1 ^> 



Ml 



> 



i - mi V 1 - y 



^ In 



Mi 



1 - Mi 



> Aln 



and so 



In 



Ai( M i,y) = <( MA) 2/ 5 

+oo if y < I 
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For 1 < A < X%: 



i 



(l-y)V i?A-l" A i?A-l V(l-y) A ' 1 " A 



Rx-1 

where d\(y, fix) — y In [Jjj^j + (1 — v) m J ■ Rearranging we can write 

d A (y,/i 1 ) = A[2/ln(y) + (l-y) ln(l-y)] - [y ln( Ml ) + (1 - y) ln(l - w )] 

which is an afhne function of A with negative slope (y ln(y) + (1 — y) ln(l — y) < for all y £ (0, 1)) and 
di(y, /ii) — K (y, fii) > 0. Hence, for fixed < y < fii < 1 this function is positive whenever 

, , yln(^i) + (l-y) ln(l-/xi) 
< yln(y) + (l-y)ln(l-y) - 

Clearly, A2 (yui , y) > 1 and we choose Ao = min(Ai, A2). After some calculation one can show that A2 < Ai, 
and therefore that 

(^2 + <5) In + (1 - M2 - 6) In 
To obtain the constants used in the statement of the lemma we define dx ^^ := d\(y,/j,i) 

Ca, M1iM2 := Ca (1 S)~ x ^ _ R - 

This concludes the proof. 



4 Experiments 

We illustrate here the performance of Thompson Sampling on numerical experiments with Bernoulli 
rewards. First we compare in terms of cumulative regret up to horizon T = 10000 Thompson Sampling 
to UCB, KL-UCB and Bayes-UCB in two different two-arms problem, one with small rewards and 
the other with high rewards, with different gaps between the parameters of the arms. Figure [l] shows 
Thompson Sampling always outperforms KL-UCB and also Bayes-UCB for large horizons. The three 
optimal policies are significantly better than UCB, even for small horizons. 

Figure [2] displays for several algorithms an estimation of the distribution of the cumulative regret 
based on N — 50000 trials, for a horizon T = 20000 in a 10-armed bandit problem with 

H = (0.1, 0.05, 0.05, 0.05, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01). 

The first two algorithms are variants of UCB. Of these the UCB-V algorithm is close to the index policy 
to which Thompson Sampling is compared in [3] in the Bernoulli setting, but this policy is not known to 
be optimal. This algorithm incorporates an estimation of the variance of the rewards in the index which 
is defined to be, for an arm that have produced k rewards in n draws, 

k + h\og{t) k 7 ~k\ | 31og(t) 
n y n n \ n J n 

The other algorithms displayed in Figure [2] have a mean regret closer (sometimes smaller) than the lower 
bound (which is only asymptotic), and among them, Thompson is the best. It is also the easiest optimal 
policy to implement, since the optimization problem solved in KL-UCB and even the computation of the 
quantiles in Bayes-UCB are more costly than producing one sample from the posterior for each arm. 
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Figure 1: Cumulated regret for the two armed-bandit problem with \x\ = 0.2, /j, 2 = 0.25 (left) and 
/ii = 0.8, /12 = 0.9 (right). Regret is estimated as an average over N — 20000 trials. 



UCB UCB-V DMED 




10 2 10 3 10 4 10 2 10 3 10 4 10 2 10 3 10 4 



KL-UCB 

500 i ■ ■ — 



400 




10 2 10 3 10 4 

n (log scale) 



Bayes-UCB 

500 i ■ 

400 
300 




10 2 10 3 10 4 

n (log scale) 



Thompson 

500 1 ■ ■ ■ 

400 




10 2 10 3 10 4 

n (log scale) 



Figure 2: Regret of the various algorithms as a function of time (on a log scale). On each graph, the red 
dashed line shows the lower bound, the solid bold curve corresponds to the mean regret while the dark 
and light shaded regions show respectively the central 99% and the upper 0.05% 
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5 Discussion 



This paper provides the first proof of the asymptotic optimality of Thompson Sampling for Bernoulli 
bandits. Moreover the proof consists in a finite time analysis comparable with that of other known optimal 
policies. We also provide here simulations showing that Thompson Sampling outperforms currently known 
optimal policies. 

Our proof of optimality borrows some ideas from Agrawal and Goyal's paper, such as the notion of 
saturated arms. However we make use of ideas together with our own to obtain a stronger result, namely 
control over the tail of Ni t t rather than its expectation. This is a valuable result which justifies the 
complexity of the proof of Proposition 2. Indeed control over these tails allows us to give a simpler finite 
time analysis for Thompson Sampling which is closer to the arguments for UCB-like algorithms, and also 
yields the optimal asymptotic rate of Lai and Robbins. 

Thanks to the generalisation pointed out in pQ, the Bernoulli version of Thompson Sampling can 
be applied to bandit problems with bounded rewards, and is therefore an excellent alternative to UCB 
policies. It would also be very natural to generalise Thompson to more complex reward distributions, 
choosing a prior appropriate for the assumptions on these distributions. Indeed, even in complex settings 
where the prior is not computable, Thompson Sampling only requires one sample from the posterior, which 
can be obtained efficiently using MCMC. Encouraging numerical experiments for reward distributions in 
the exponential family using a conjugate prior suggest that a generalisation of the proof is achievable. 
However this poses quite a challenge since the proof here is often heavily dependent on specific properties 
of Beta distributions. A natural generalisation would need a prior-dependent finite-time result controlling 
the tail probabilities of posterior distributions as the number of samples increases. 
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